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^^V\'e investigate laser induced nonsequential multiple ionization using a simple statistical model in 
j^which an electron recoUides inelastically with its parent ion. In this collision, it thermalizes with 

«the remaining A'^ — 1 bound electrons within a time interval At. Subsequently, all the A'^ electrons 
'^eave. We address the question of how the above time delays influence the individual contributions 
^Ijfrom the orbits in which the first electron, upon return, is accelerated or decelerated by the field, 
(^Respectively, to the ion momentum distributions. In both cases, the time delays modify the drift 
^icnomenta obtained by the electrons when they reach the continuum at t + At, by moving such 
(—( times towards or away from a crossing of the electric field. The contributions from both types of 
J— j;ollisions are influenced in distinct ways, and the interplay between such trends determines the widths 
l_^nd the peak momenta of the distributions. Specifically in the few-cycle pulse case, we also show 

that such time delays do not affect the shapes of the momentum distributions in a radical fashion. 

pence, even with a thermalization time At, nonsequential multiple ionization could in principle be 
^used for absolute-phase diagnosis. 

(N ■ 
^ ■ 

^""1 Introduction 
l> ■ 

v^The physics of non-sequential double ionization in rare gases by low-frequency, 
Cintense laser fields is, to a very large extent, understood. In particular for neon, 
c/3the electron momentum distributions with peaks at non-vanishing momentum 
.Components pi\\ = P2\\ = -iI.2^JlJp, along the laser-field polarization, or the 
^on momentum distributions peaked at the parallel momenta P\\ = ib4y^f^, 
cohere Up denotes the ponderomotive energy, can be explained by a laser- 
Gnduced inelastic recollision process [1-3]. Specifically, the first electron leaves 
J>an atom by tunneling ionization at a time t', propagates in the continuum, 
Iljpbeing accelerated by the field and, subsequently, at a time t, recollides with its 
*''^^arent ion. The second electron is then freed by electron-impact ionization [4]. 
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This process has been successfully dealt with using an S-Matrix approach 
and its classical counterpart [5-9]. In particular, if the laser-field intensity is 
high enough, both quantum-mechanical and classical models yield very simi- 
lar momentum distributions"^. Wc have shown that this agreement occurs for 
several types of electron-electron interaction, and in the absence or presence 
of Coulomb repulsion in the electron final states [6-8]. Furthermore, we have 
verified that the outcome of the classical and quantum mechanical computa- 
tions also exhibit a striking agreement if the driving field is a few-cycle laser 
pulse [9, 12]. In particular, both the classical and quantum mechanical distri- 
butions are strongly dependent on the so-called absolute phase, i.e., the phase 
difference between the pulse envelope and its carrier oscillation. We have used 
this fact in order to propose a scheme for diagnosing the absolute phase [12], 
which has also been experimentally realized [13]. 

If one is dealing, however, with more than two active electrons, there are 
many open questions concerning the physics behind laser-induced multiple 
ionization. This is in particular true for triple and quadruple ionization of 
neon by near infra-red (a; = 0.057 a.u.) laser fields in an intensity range be- 
low 2 X 10"'^''' W/cm^ [2, 14, 15]. The distributions observed in such experi- 
ments, approximately peaked at non-vanishing ion momenta = ±2Ny/Tjp, 
where N denotes the number of electrons involved, hint at a non-sequential 
physical mechanism. The theoretical modeling of non-sequential multi-electron 
processes for N > 2, however, poses a far greater challenge than in the two- 
electron case. For instance, already for three electrons, which is the simplest 
scenario to be taken into consideration, one must include at least six Feyn- 
man diagrams when computing the corresponding transition amplitude, each 
of which is rather cumbersome to implement. 

In order to tackle such a problem, in previous work [16], we constructed 
a simple thermalization model, similar to those employed in, for instance, 
nuclear [17] or molecular [18] physics. We assumed that the first electron, 
upon return, recollides with its parent ion and shares its kinetic energy with 
N — 1 bound electrons. All the N electrons are then freed after a time interval 
At subsequent to the recollision time t, which, to first approximation, is taken 
to be constant. We regard this time as the upper bound for a thermalization 
time, which is necessary for the kinetic energy to be redistributed among the 
N electrons. 

The above-mentioned process takes place within a fraction of the cycle of the 
laser field. Since the duration of a typical titanium-sapphire pulse, employed 
in such experiments, is of the order of 2.7 fs, At falls within the attosecond 



This, however, no longer holds if tlie intensity is near threshold, i.e., just enough for the second 
electron to be released. For a detailed discussion c.f. [10], and for measurements of nonsequential 
double ionization below the threshold see, e.g., [11]. 
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regime. Indeed, recently, by comparing our thermalization model with the 
available experimental data for triple and fourfold ionization of neon [14,15], we 
have been able to estimate thermalization times of the order of 460 attoseconds 
[16] . Using such a time delay, wc computed ion momentum distributions whose 
widths and peak momenta agreed very well with those in the experiments. 

To a very large extent, our thermalization model contains classical ingre- 
dients and is a generalization of the classical model already employed in the 
two-electron case [7,8]. In [16], we have also observed that the peaks and the 
width of the momentum distributions depend very strongly on the time de- 
lay At. The position of the peaks is a direct consequence of the kinematics 
involved in the process, as the final drift momentum the N electrons acquire 
from the field, will strongly depend on the time t + At when they reach the 
continuum. 

The above-stated physical mechanism, however, is not completely agreed 
upon. Indeed, it could well be that a broadening in the electron momentum 
distributions and the displacement in the peak momenta are caused by, for 
instance, partially sequential ionization processes [15], or excitation tunnel- 
ing mechanisms [19] , in which, during the rccoUision, the electrons undergo 
a transition to an excited state, from which they subsequently tunnel (for 
an example of such mechanisms for different atomic species see, e.g., [20]). 
Specifically for large thermalization times, the distributions are very much 
concentrated near vanishing parallel momenta, and resemble those obtained 
with excitation-tunneling. Apart from that, even assuming that thermaliza- 
tion through electron-impact ionization is the only or at least the dominant 
physical mechanism, it is not clear how the time delay At affects the width of 
the distributions. 

Another open question concerns the infiuence of the types of collisions on 
the momentum distributions. In fact, for the recombination or rescattering 
phenomena occurring in the context of atoms in strong laser fields, there exist 
two main orbits along which an electron may return to its parent ion with the 
same kinetic energy. Such orbits are known as "the long orbit" , or "the short 
orbit", and have been extensively investigated in the context of high-order 
harmonic generation and above-threshold ionization [21] and, more recently, 
for nonsequential double ionization [22, 23]. An electron following the long 
orbit leaves the atom at a time t' < O.ST, where T = 2tt/uj denotes a field 
cycle, and returns after the minimum of the electric field, whereas an electron 
along the short orbit reaches the continuum at t' > O.ST and returns before 
the minimum of the electric field. Hence, in the former and in the latter case 
the returning electron is decelerated and accelerated by the field, respectively. 

In recent years, both orbits have been investigated in the context of non- 
sequential double ionization using classical ensemble models [22]. In particular, 
they have been associated to two main types of collisions. The orbits for which 
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the electron is being decelerated have been related to the so-called "slow-down 
collisions" , whereas the short orbit has been linked to the so-called "speed-up 
collisions". Depending on several parameters, such as the field intensity or 
frequency, one type of collisions may provide the dominant contributions to 
the yield, or both types may be equally important [27]. For instance, for the 
parameters in [22], it was shown that the slow-down collisions dominate. 

In this paper, we address the question of how the contributions from the 
short and the long orbit to the multiple-ionization momentum distributions 
are affected by the time delay At, for triple and fourfold ionization. We are 
mainly concerned with the influence of both sets of orbits on the width and 
the peaks of the ion momentum distributions, for resolved and non-resolved 
ion momentum components transverse to the laser-field polarization. We take 
the driving field to be a monochromatic wave or a few-cycle pulse, and make 
a detailed analysis of the similarities and differences between both cases. This 
manuscript is organized as follows: In the subsequent section, we provide a 
brief discussion of our thermalization model, with the explicit expressions for 
the ion momentum distributions. In Sec. El we display our results for the 
different sets of orbits, together with the total yields and, finally, in Sec. 01 we 
state our conclusions. 

2 Model 

We will briefly recall the statistical model employed in our previous publica- 
tion [16], in order to determine the momentum distributions in N-fold non- 
sequential ionization. This model is mainly composed by three main ingredi- 
ents. 

Firstly, we assume that one electron tunnels into the continuum at an instant 
t' according to the quasi-static rate [25]. 



where E{t') and \Eqi \ denote the electric field at the time the electron left and 
the first ionization potential, respectively. 

Subsequently, this electron propagates under the sole influence of the laser 
field, and obeys the classical equations of motion. At a later time t, it may be 
driven back towards its parent ion, recolliding inelastically with it. The return 
condition ri(t) = for the first electron yields 



Rit') ~ |^(t')r'exp -2i2\Eoi\fyi3\E{t')\) 



(1) 




(2) 



where A(t) is the vector potential, and F{t) = J* A{T)dT. Eq. ((JJ provides a 
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clear geometrical picture for the return condition: for a given t', the return 
time t is given by the intersection of F{t) with its tangent at t' [26]. We will 
employ this construction in Sec. 01 when discussing our results (c.f. Fig. 2). 

Upon recollision, the energy E^^t{t) = [A{t) - A{t')f/2 of the returning 
electron is completely thermalized among such an electron and the — 1 
electrons to be freed. After the time interval At, the distribution of energy 
and momentum over the electrons is only governed by the available phase 
space. We will consider here that the time delay At is constant. Finally, at the 
time t + At, the electrons reach the continuum with the total kinetic energy 
E'ret — The quantity e'^'' = X]^=2l-^0n| denotes the total ionization 

potential of the — 1 (up to the recollision time t inactive) electrons. 

The distribution of final electron momenta p„ (n = 1, . . . , A^) is proportional 

to 

F(pi,p2, . . . ,P7v) = j dt'R{t')6 (eJ^^) - E,,t{t) 

1 ^ \ 
+-^[p„ + A(t + At)f , (3) 

n=l ) 



where p„, and A(t + At) denote the final electron momenta and the vector 
potential at the time the electron leaves, respectively. The integral extends 
over the ionization time t'. The b function expresses the fact that the total 
kinetic energy of the A'^ participating electrons is fixed by the first-ionized 
electron at its recollision time t. 

In the 3 — A^ dimensional momentum space (pi, . . . , Pat), the argument of 
the b function corresponds to the equation of a hypersphere 



- Y^PnW + A(t + At)]2 + - ^ p^^ = Sret(t) " Ef\ (4) 
n=l n=l 



where and Pn±(n. = 1, ■■■,N) denote the electron momentum components 
parallel and perpendicular to the laser-field polarization, centered at —A{t + 
At), and whose radius is determined by the difference between the kinetic 

energy E^ct{t) of the first electron upon return and the absolute value E^^^ of 
the total binding energy which must be overcome. Within the region delimited 
by this hypersphere, in a classical framework, the process we are dealing with 
takes place. Furthermore, if the transverse momenta are kept constant, Eq. 
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(jlj) can be written as 



- Y,[Pn\\ + A{t + At)f = Erct(t) - El 



(5) 



n=l 



where Eq = Eq + ^ Z]^=i Pn± '^^'^ regarded as an effective binding 
energy wfiicli the N electrons must overcome. Such an energy depends on 
the transverse momentum components, and is minimal when they vanish. By 
employing an adequate choice of field and atomic parameters, such as, for 
instance, a few-cycle pulse, one can manipulate the radius in making 
whole momentum regions collapse or appear. In [9,12], we have exploited this 
in the two-electron case in order to determine the absolute phase. 

The only free parameter of this model is the time delay At between the rec- 
ollision time and the time when multiple ionization occurs. It is the sum of the 
time it takes to establish the statistical ensemble, i.e., the thermalization time, 
and a possible additional "dwell time", until the electrons become free. This 
model is an extension to NSMI of a classical model introduced for NSDI in 
Refs. [7,8] for At = 0. Sufficiently high above threshold, it produced momen- 
tum distributions that were virtually indistinguishable from their quantum- 
mechanical counterparts. 

In order to integrate over unobserved momentum components, we exponen- 
tialize the 6 function in Eq. ^ using its Fourier representation 



Infinite integrations over the momenta p„ can then be done by Gaussian 
quadrature. The remaining integration over the variable A is performed em- 




(6) 



ploying [28] 





where x'^ = x'^9{x), with 9{x) the unit step function and e — > +0. 
The momentum distribution of the ion is then 



F(P)= / [JdW P + 5^Pn p(pi,P2,...,Piv) 




, 3N 3 

(27r) — "2 



/ 



3fl 

dt'R{t') (ASjv,ion)+^ 



2 



(8) 



Af3/2r(3(iV - l)/2) 
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with A^AT^ion = E,et{t) - E^^^^ - ^[P - NA{t + At)f, where P gives the 
momentum of the ion. In ©, we have used the momentum conservation P = 
— '^n=i Pn, which is vahd in case the photon momenta can be neglected. 

If the ion momentum component perpendicular to the laser polarization 
is entirely integrated over, the remaining distribution of the longitudinal ion 
momentum Pii reads 



iVr((3iV - l)/2) 



I dt'R{t'){AEN,ion\\)J K (9) 



where AF;v,ion|| = ^ret(t) - E^/"^ + Ai)]^ 

If just one transverse-momentum component (for instance, -Pi, 2) is inte- 
grated while the other one (-Pi,i = -Pi) is observed, the corresponding distri- 
bution is 



\ —1 

^^""^ ' ''dt'R{t'){AE^^,,,^^^)f-' (10) 



NT{3N/2 - 1) 



with A-Ejv,ion|l± = A-Bion|| — ^P]_- In the following section we will explicitly 
compute ion momentum distributions employing our thermalization model, 
i.e., from Eqs. Q and (fTTH) . 



3 Momentum distributions 
3.1 Monochromatic fields 

As a starting point, we will consider a monochromatic laser field E(i) = 
—dA(t)/dt, for which the vector potential is of the form 



A(t) = ^ocos(u;t)e2, (11) 

with frequency uj, amplitude Aq and polarization axis €z . This is a good ap- 
proximation for long pulses. In this case, A-E'^ion||(i-||, t) = Ai^'^iojj||(— P||,^-|- 
^/2), where T = 2tt/uj denotes a period of the driving field, so that the mo- 
mentum distributions are symmetric with respect to P|| = 0. Clearly, this also 
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Figure 1. (Color online) Distribution of the longitudinal ion momentum for triple nonsequential 
ionization of neon at 7 = 2 X lO^^W/cm^, for time delays At = 0, At = 0.05T, At = 0.17T, and 
At = 0.25T, together with individual contributions from the long and the short orbits. The ion 
momenta are given in units of [Up]^^^ 



holds for A£'^ion||^. In our computations, we consider the shortest pair of 
trajectories along which the first electron returns to its parent ion, and, in 
particular, the individual contributions from the longer and the shorter orbit 
within such a pair. Such orbits coalesce at the boundary of the region deter- 
mined by the hypershpere (jU for which, classically, nonsequential multiple 
ionization is allowed to occur. In the particular model employed in this paper, 
such a boundary is determined by the minimal momenta the electrons need 
in order to overcome the binding energy E^^'^ , and the maximal momenta for 
which the first electron is able to return. 

In Fig. ^ we display the yield computed assuming that the first electron 
returns to its parent ion either along the long or the short orbit, together 
with the total momentum distributions, for triple ionization of neon. In the 
upper panel [Fig. ^(a)], the thermalization time is taken to be vanishing. As 
an overall feature, the total yield is approximately peaked at = ib6y^Up, 
as expected from classical arguments. Furthermore, the shorter orbit yields 
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distributions which are shghtly shifted towards lower momenta, in comparison 
to those from the longer orbit. Such a shift is corrected if a time delay At = 
0.05T is introduced [Fig. ^(b)]. This can be easily understood in terms of 
the momentum transfer from the field to the electrons. Such a momentum 
transfer is maximal if all electrons leave at a crossing of the field, i.e., for 
t + /S.t = 71 vr. For vanishing time delay, the electrons released by a speed up 
collision, i.e., by an electron returning along the short orbit, leave at roughly 
ts = 0.95T, while for the slow-down collisions, this happens at = LOST. 
The time delay At = 0.05T moves the release times tg + At for the short-orbit 
contributions towards the crossing and their long-orbit counterparts tg + At 
away from it. As a direct consequence, the contributions from the former set 
of orbits move towards higher absolute momenta, while those from the latter 
set move towards lower absolute momenta. 

If the time delay increases [Fig. ^(c)], the peak momenta from the short- 
orbit contributions starts to move again towards lower absolute values. This is 
due to the fact that, for At > O.OST, the time ts + At all electrons leave starts 
to move away from the crossing of the driving field. Hence, the drift momenta 
of the electrons freed by the speed-up collisions decrease. One should note, 
however, that the times t^ -|- At, for the electrons released by the slow-down 
collisions, also move farther away from the crossing. In this latter case, how- 
ever, the peak momenta decrease in a far more radical way. The fact that the 
two individual contributions are now peaked at momenta which are very far 
apart leads to a broadening in the momentum distributions. In this context, 
it is worth mentioning that the time delay At = O.ITT yields the best agree- 
ment with the experimental findings, for nonsequential triple and quadruple 
ionization [14,15]. This issue is discussed in more detail in [16]. 

For exceptionally large time delays, such as those depicted in Fig. ^(d), 
the contributions from the short orbit move back to lower momenta, and those 
from the long orbit overlap at vanishing ion momenta. As a direct consequence, 
the total yield merges in a single peak at vanishing momentum. Such features 
resemble to a very large extent those obtained if one considers sequential 
physical mechanisms, or the situation where the electrons undergo a transition 
to an excited state, and, subsequently, tunnel into the continuum. One should 
note, however, that for such a parameter range, there is poor agreement with 
the experimental results. 

The effect observed in Fig. ^ is shown in more detail in Fig. |^ using the 
tangent construction mentioned in the previous section. In the figure, we con- 
sider the return condition ^ for the short and long orbits. For vanishing time 
delays, the electron starting at t'^ and t^ releases the remaining electrons at 
the times tg before or t^, after the crossing of the field, respectively. Such times 
are determined by the intersection of the tangents Q with the dashed line. If, 
however, there is a time delay At, the time t + At in which all electrons leave 
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Figure 2. (Color online) Start and return times for the long and short orbits, denoted by 
(*Siis)i respectively, for vanishing and non- vanishing time delays At. 

will be determined by the intersection of the tangents © with the solid line. 

In Fig. |51 we depict the contributions from the slow-down and the speed-up 
collisions, as well as from the total yield, for the quadruple-ionization case and 
the same time delays as in Fig. ^ They exhibit the same overall behavior as 
in Fig. ^ the peak momenta for the short-orbit contributions moves towards 
higher absolute values for At = 0.05T [Fig. 01 (b)], and, as At is further 
increased, the peak momenta start to decrease [Fig. 01 (c)]. The maxima of 
the long-orbit contributions always shift towards Py = 0, for increasing At, 
as shown in Fig. ini(b)-(d). This behavior is expected, since within our model, 
both triple and quadruple nonsequential ionization are governed by the same 
physical mechanism, namely energy transfer to all active electrons due to the 
recoUision of the first electron with its parent ion. The main difference is 
the fact that all electrons must overcome a larger energy in order to 
reach the continuum. Therefore, for the same laser- field intensity, as the charge 
state increases, the radius of the hypersphere Q decreases and so does the 
momentum region for which the process in question occurs. Furthermore, the 
peak momenta should increase to the vicinity of P|| = itSy^L^. This leads, in 
general, to distributions which are much more localized around the maxima 
than in the triple-ionization case. 

For a time delay At = 0.25T, however, instead of a single peak at vanishing 
momenta, the total distribution exhibits a plateau, with humps near its low- 
and high-energy ends [Fig.|31.(d)] . This is due to the fact that, while the region 
near P^^ = is filled by the contributions from the long orbit, the contributions 
from the short orbit are still at a much larger momenta than in the two electron 
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Figure 3. (Color online) Distribution of the longitudinal ion momentum for quadruple 
nonsequential ionization of neon at / = 2 X lO^^W/cm^, for time delays A.t = 0, At = 0.05T, 
At = 0.17T, and At = 0.25T, together with individual contributions from the long and the short 
orbits. The ion momenta are given in units of [C/p]^/^. 

case, and still cause such humps. However, the distributions from the negative 
and positive momentum regions still merge, as in the triple-ionization case. 

Finally, in Fig. |1J we analyze the individual contributions of both orbits if 
the transverse momenta are only partly integrated over, for time delays Ai = 
and Ai = O.ITT. The distributions from the short and the long orbits exhibit 
a displacement in their peak momenta which is similar to those obtained in 
the non-resolved case. In particular, the elongation observed in experiments 
for such distributions is directly related to this displacement. 



3.2 Few-cycle pulses 

In this section, we will take the driving field to be a few-cycle pulse E(f) = 
—dA{t)/dt, corresponding to the vector potential 



A{t) = Aosm'^[ujt/{2n)]sm{ujt + (t>)ez 



(12) 
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Figure 4. Ion momentum distribution of nonsequential triple ionization of neon at 
/ = 1.5 X lO-'^^W/cm^ , and time delay At = O.ITT. Panels (a), (b) and (c) correspond to the 
contributions of the short orbit, of the long orbit, and to the total yield, respectively. In panel (d), 
the total yield for vanishing time delay is given for comparison. The ion momenta are given in 

atomic units. 



where 4> is the absolute phase, and n is the number of cycles. We choose a 
four-cycle pulse (n = 4), and consider triple and quadruple ionization. Apart 
from their wide range of applications (for reviews see, e.g., [24]), few-cycle 
pulses are very attractive in the context of nonsequential multiple ionization. 
This is due to the fact that several physical mechanisms which would compete 
with electron-impact ionization, such as excitation-tunneling ionization, are 
strongly suppressed in the few-cycle case [14]. Therefore, in principle, one 
expects cleaner measurements. 

In Fig. El we present ion momentum distributions for triple ionization of 
neon, computed with the driving field (|12j) . As a starting point, we take the 
thermalization time At to be vanishing. In contrast to the monochromatic 
case, such distributions are in general asymmetric with respect to Py = 0. 
This is expected since, for a few cycle pulse, A{t) = ±A{t + nT/2) holds, 
if at all, only for very specific sets of parameters. Therefore the relation 
Ai?jv,ion|| ) = ^^N,ion\\i~P\\j't + T/2) IS in general not fulfilled. Depend- 
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Figure 5. Ion momentum distribution of nonsequential triple ionization of neon by a four cycle 
pulse of intensity / = 2 X IQ-'^^W/cm'^, and several absolute phases </>. We assume that the electrons 
leave instantaneously after recoUision {At = 0). The ion momenta are given in units of [f/p]^/^. 



ing on the absolute phase, the distributions are strongly localized in either the 
positive or negative momentum region. Around a critical phase, they start to 
shift from one momentum region to the other. 

Such a behavior is very similar to that observed for nonsequential double 
ionization (N = 2). In the two-electron case, this shift has been employed 
to diagnose the absolute phase [9, 12, 13]. Furthermore, it has been explained 
as a change in the dominant set of trajectories for an electron rescattering 
with its parent ion. Such a set is determined by the interplay between the 
rate with which the first electron is released in the continuum, and the 
available phase space: if there was a large probability for the first electron to 
leave along a particular set of orbits, and if it returned with enough energy 
to cause appreciable electron-impact ionization, one expects the contributions 
from such orbits to the yield to be large. 

Within the present model, this explanation also holds for N > 2, and is 
illustrated in Fig. [S] for the parameter set of Fig. |SJ In the upper and lower 
panels, we present the quasi-static tunneling rate and parallel momentum com- 
ponents for the residual ion, respectively, as functions of the emission time t'. 
In the latter case, we consider the transverse ionic momentum components to 
be vanishing. This gives an upper bound for the momentum region for which 
electron-impact ionization is allowed to occur, since, in Eq. ©, the effective 



14 Speed-up and slow-down collisions in laser-induced nonsequential multiple ionization 



binding energy Eq is minimal. 

In the figure, we identify at least three sets of orbits, whose start times t' 
are highly dependent on the carrier-envelope phase 4>. For instance, for (p = 0, 
there exist three such pairs: (1,2), (3,4) and (5,6). Thereby, the odd and 
even numbers refer to the longer and the shorter orbit in a pair, respectively. 
Orbits (1,2) and (5,6) contribute to the yield in the positive momentum re- 
gion, whereas (3,4) leads to a peak in the region Py < 0. Since, for (1,2), the 
quasi-static rate is very small, the contributions from this set of orbits is 
negligible. Hence, the features observed in Fig. 5 will be determined by the 
remaining two pairs. In fact, orbits (3, 4) and (5, 6) lead to the large peak for 
negative momenta and to the small peak for positive momenta, respectively. 
For the set (3, 4), there is a large momentum region for which the first electron 
may dislodge the remaining two, and a non-negligible ionization rate. There- 
fore, its contributions to the yield are expected to be quite prominent. For the 
set (5,6), even though the tunneling rate is very large, there is only a small 
region for which electron-impact ionization is allowed. Hence, the pair (3, 4) 
dominates. 

This picture starts to change around the critical phase (pc = 0.37r, for which 
the yield is approximately symmetric. In comparison to the case 4> = 0, there 
is now a decrease in the quasi-static rate for (3, 4) and an increase in the mo- 
mentum region for (5,6). Therefore, the contributions from the former orbits 
decrease and those from the latter orbits increase, respectively. Beyond the 
critical phase (for instance, for (p = O.Svr), the yield is mainly concentrated in 
the positive momentum region, since now (5, 6) prevail. For this latter case, 
there is also a further set of orbits which contribute to the yield in the region 
P|| < 0, namely (7,8). For such a set, the rate is very large. However, the 
region for which electron-impact ionization is allowed to occur is very small 
and (5, 6) is still the dominant pair of orbits. 

For a higher charge state, the very same line of argumentation applies. The 
main difference is that the first electron will need to overcome a larger binding 
energy in order to release the remaining electrons. As a direct consequence, the 
radius of the hypersphere (jlj) will be smaller so that whole momentum regions 
may collapse, or contribute in a much less pronounced way to the yield. For 
instance, for nonsequential quadruple ionization, for the absolute phase (p = 
the small peak on the positive momentum region is absent. This is due to the 
fact that the momentum space for the orbits (5, 6) is much smaller than in 
the three-electron case, so that their contributions to the yield are negligible. 
For (p = O.Svr, we have also observed that the small peak in the negative 
momentum region is not present, for the same reason. 

We now consider non-vanishing thermalization times and a constant abso- 
lute phase, which has been chosen so that the yield is almost entirely concen- 
trated in the negative momentum region. The results of these computations 
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Figure 6. (Color online) Quasi-static tunneling rate (upper panels), and parallel momenta Py 
(lower panels), as functions of the start times t', for triple and quadruple ionization of neon. We 
consider absolute phases <f) = 0, (p = O.Stt and (f> = O.Stt, the same driving-field parameters as in the 
previous figure, vanishing transverse momenta P± and thermalization times At = 0. The numbers 
(i, j) in the figure indicate a pair of orbits, with the odd and even numbers referring to the longer 
and the shorter orbit in a pair, respectively. The start time t' is given in units of the field cycle and 

the ion momenta in units of [f/p]^/^. 

are displayed in Fig. d As an overall feature, the time delays At influence the 
width and the position of the peaks of the distributions. They do not modify, 
however, their asymmetric shapes or the momentum region the yield is con- 
centrated. This is due to the fact that a non-vanishing At changes the drift 
momenta the N electrons acquire from the external laser field when they reach 
the continuum, i.e., the left-hand side in Eqs. (@J) and ©, but does not affect 
the quasi-static tunneling rate or the binding energy the electrons must 
overcome [i.e., the right-hand side in Eqs. (@J, ©]• Therefore the dominant 
pairs of orbits remain the same. In other words, by changing At one influences 
the center of the hypersphere (jU, but not its radius. 

We shall now have a closer look at how the width of the distributions is 
affected by the interplay between the long and the short orbit. In Fig. |H1 
we display the ion momentum distributions for quadruple ionization in the 
few-cycle pulse case, together with the individual contributions from the long 
and the short orbits. For a few-cycle pulse, we define such orbits with respect 
to the pair which gives dominant contributions to the yield. For this specific 
parameter set, this means orbits (3, 4) in Fig. |H1 The long- and short-orbit 
contributions, i.e., the individual contributions from orbits 3 and 4, respec- 
tively, behave in a similar way to the monochromatic-field case (Figs. ^ and 
131) • Indeed, for vanishing thermalization times, the momentum distributions 
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Figure 7. (Color online) Distribution of the longitudinal ion momentum for triple [panel (a)] and 
quadruple [panel (b)] nonsequential ionization of neon, for a four-cycle pulse of intensity at 
I = 2 X lO^^W/cm^, absolute phase </> = 0, and several time delays. The momenta are given in 

units of [r/p]i/2. 



from the slow-down collisions are peaked at a slightly larger absolute momen- 
tum than those from the short orbit. For At = O.IT [Fig. IHl(b)], the former 
and the latter contributions have moved away from or towards Py = 0, re- 
spectively. Once more, as the thermalization time increases, the peaks of both 
contributions shift towards vanishing ion momenta, as shown in Fig. |Hl(c). 

A qualitative difference, however, between Fig. |H1 and its monochromatic 
counterpart is that, for large time delays, one no longer obtains a plateau 
in the vicinity of vanishing ion momenta, but a single peak in the negative 
momentum region. This is due to the fact that, for a monochromatic driving 
field, the symmetric peaks in the positive and negative momentum regions 
start to merge. For a few-cycle driving pulse, on the other hand, by taking an 
appropriate absolute phase, the distributions can be chosen in such a way that 
they are concentrated only in one momentum region. Hence, it is possible to 
avoid that a plateau or a single peak occurs near P\\ = 0. 



4 Conclusions 



We have investigated the influence of the shorter and the longer orbit of an 
electron recolliding inelastically with its parent ion on ion momentum distribu- 
tions, within the context of a thermalization model for laser-induced A^— fold 
nonsequential ionization. In this model, the first electron, upon recollision, 
shares the kinetic energy it acquired from the external laser field with the 
remaining — 1 bound electrons. All A^ electrons then leave after a time in- 
terval At. We have taken the driving field to be a monochromatic wave and a 
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Figure 8. (Color online) Distribution of the longitudinal ion momentum quadruple nonsequential 
ionization of neon, together with the individual contributions from the long and short orbits, for a 
four-cycle pulse of intensity / = 2 X lO^^W/cm'^, absolute phase (j> = 0. We consider the time delays 
M = 0, At = O.IT and At = 0.2T. The ion momenta are given in units of [f/p]^/^. 



few-cycle laser pulse. 

We have traced the variations in the widths and peaks of such momentum 
distributions, with respect to the time delays At, back to the momentum 
transfer from the field to an electron traveling along the short or the long orbit, 
and recolliding with its parent ion. The momentum transfer is maximal at a 
crossing of the electric field. With increasing time delays, the time for which 
the electrons released by a slow-down collision, i.e., by an electron moving 
along the long orbit, will leave their parent ion moves away from this crossing. 
As a direct consequence, the momentum transfer will decrease and the peak 
momenta will move towards lower absolute values. The situation is however 
quite different for the electrons released by a speed up collision. In this case, 
following the increase of time delay, the release time t + At will initially (i.e., 
for At < 0.05T) approach a crossing, and subsequently (i.e., for At > 0.05T) 
move away from it. This will make the peaks of the pertinent momentum 
distributions move to larger absolute values, and then back towards = 0. 
The width of the resulting total momentum distributions, as well as their 
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peaks, will result from the interplay between such trends. This also explains 
the elongated patterns observed for the momentum distributions, in case the 
perpendicular momentum of the ion is resolved. 

In both monochromatic and few-cycle cases, we considered the dominant 
pair of orbits, which, for each situation, is chosen in a slightly different way. 
For monochromatic fields, we took the two shortest trajectories, i.e., those for 
which the first electron spends the shortest time in the continuum. Physically, 
this pair is dominant, since the spreading of the electronic wave packet is less 
pronounced than for longer pairs of orbits. Hence, the overlap with bound- 
state wave packets will be appreciable^. In the few-cycle case, the dominant 
pair is determined by the interplay between the tunneling rate for the 
first electron and the available momentum space. A pair of orbits will only 
contribute in a relevant way to the yield if the first electron is ejected with 
a large probability and if, upon return, it has enough energy to release the 
second electron (see, e.g. [9,12], for more details). 

Specifically for few-cycle pulses, we have shown that non- vanishing thermal- 
ization times do not play a very important role in determining the shape of 
the momentum distributions. This is due to the fact that a time delay At 7^ 
between recollision and multiple ionization only alters the drift momenta the 
N electrons acquire from the field when they reach the continuum. However, 
this delay does not modify the key ingredients which determine the dominant 
pairs of orbits and therefore the shape of the momentum distributions: the 
tunneling rate for the first electron and its kinetic energy upon return. In par- 
ticular, both ingredients are highly dependent on the carrier-envelope phase, 
so that the shape of the momentum distributions can be used for measuring 
this parameter [12,13]. This means that, even in case there is a non- vanishing 
thermalization time, one could, in principle, still employ nonsequential mul- 
tiple ionization for absolute-phase diagnosis. Finally, few-cycle driving pulses 
may be very useful for testing and validating the thermalization model fur- 
ther, since many other ionization channels, such as excitation-tunneling or 
sequential ionization, are strongly suppressed in this case. 
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Rudenko for discussions on nonsequential multiple ionization with few-cycle 
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^We would like to stress out that wave-packet spreading is not included in our classical model. It 
is, however, present in the quantum-mechanical model we have employed in the two-electron case, 
which is its quantum-mechanical counterpart. In this latter case and as expected from our physical 
intuition, we verified that the main contributions to the yield come from the shortest pair. 
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